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INTRODUCTION  : 

The  term  "consumer  durables"describes  consumer  goods  that  have  the 
characteristic  of  providing  a  constant  stream  of  services  over  the 
period  of  their  useful  life.  Examples  of  such  goods  are  household 
appliances  such  as  washing  machines  and  refrigerators.  In  goods  of  this 
type  questions  of  reliability,  durability  and  maintainability  of  the 
product  become  important,  especially  where  the  product  is  long-lived 
and  not  subject  to  rapid  obsolescence  or  vogue  change. 

The  properties  of  reliability,  durability  etc.  enter  into  the  costs 
incurred  in  providing  continuous  service  over  the  life  of  the  product. 
These  costs  are  due  to  repair,  maintenance,  servicing  and  part  replace- 
ment. Clearly  the  greater  the  reliability  of  the  product  the  less  these 
costs  will  be.  However  the  improvement  of  these  characteristics  can 
only  be  bought  at  a  price  -  that  of  increased  design  and  production 
costs  -  which  is  going  to  be  reflected  in  the  price  of  the  product 
to  the  consumer.  From  the  buyer's  point  of  view  there  exists  therefore, 
some  optimal  level  of  reliability  and  maintainability  depending  on 
his  relative  preferences  for  present  and  future  costs. 

As  for  the  manufacturer  ,  he  is  expected  to  provide  a  product  that  delivers 
a  given  stream  of  services  at  a  level  of  reliability  commensurate  with 
the  sale  price.  However,  given  that  consumers  need  not  be  accurately 
informed  about  what  future  costs  to  expect,  there  is  no  incentive  for 
manufacturers  to  be  maximally  productive  with  regard  to  product  reliability. 
That  is  to  say,  if  we  plot  on  a  graph  points  representing  technologically 
feasible  combinations  of  price  versus  future  repair  costs,  it  may  be  possible 
for  manufacturers  to  operate  at  points  far  away  from  the  production- 
possibility  frontier  which  is  formed  by  the  southwest  boundary  of  the  set 
of  feasible  combinations.  (Figure  1) 
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The  economics  of  the  situation,  and  the  behavior  of  consumers  and 
manufacturers  have  been  analysed  in  detail  by  Hausman[l]  and 
Greenwood  [2].  The  interaction  with  prices  in  markets  for  used  durable 
goods  has  been  extensively  discussed  by  Dana  [3].  The  general  area  of 
repair  and  maintainance  in  the  household  appliances  industry  has  been 
reviewed  in  [4]  together  with  a  discussion  of  significant  problem  areas 
and  policy  recommendations.  Suffice  it  to  say  that  it  is  possible  for 
diseconomies  to  exist  in  the  present  system  where  responsibility  for 
repair  and  maintainance  of  the  product  lies  entirely  with  the  consumer. 

SERVICE  CONTRACTS  .... 

Hausman  [1]  and  Greenwood  [2]  have  suggested  the  mandatory  offering  of 
service  contracts  by  manufacturers  as  an  elegant  solution  to  the  problems 
of  productivity  in  the  area  of  repair  and  maintainance  of  durables.  A 
service  contract  is  an  agreement  between  the  manufacturer  and  the  buyer 
by  which  the  manufacturer  undertakes  to  repair  and  maintain  a  product 
for  a  fee  payment  by  the  owner.  Such  an  offering  would  leave  the  customer 
the  option  of  buying  the  contract  or  choosing  to  bear  the  costs  of  repair 
etc.  himself.  In  either  case  the  result  would  be  that  the  consumer  is 
fully  informed  about  what  magnitude  of  repair  costs  he  might  expect  to 
incur.  Furthermore  the  system  benefits  from  such  an  arrangement  since 
with  the  centralization  of  the  repair  operation  in  a  factory  environment, 
economies  of  scale  can  be  realized  due  to  improved  planning  and  effective 
resource  allocation  and  utilization.  In  short  ,  a  very  good  case  can  be 
made  for  mandatory  service  contracts  from  the  viewpoint  of  consumer  pro- 
tection and  overall  system  economy. 

...  and  the  MANUFACTURER  : 


For  the  manufacturer  however,  there  is  no  distinct  advantage  to  be  gained 
from  such  an  offering  unless  he  feels  that  he  has  a  competitive  technolo- 
gical advantage  in  the  area  of  product  reliability,  that  mandatory  service 
contracts  would  allow  him  to  exploit.  Indeed  several  manufacturers  do  in 
fact  emphasise  reliability  in  their  marketing  approach  (Maytag,  VW  for 
example)  and  some  already  offer  a  service  contract  option. 


Given  that  a  manufacturer  has  decided  to  offer  a  service  contract,  it 
is  in  his  interest  to  ensure  that  he  derives  due  benefit  from  the  offer, 
since  he  is  contracting  to  provide  services  that  would  otherwise  be  provided 
by  some  other  profit  making  organization.  The  contract  should  be  seen  as 
a  product  to  be  appropriately  priced  and  marketed  and  for  which  production 
and  inventory  planning  decisions  have  to  be  made.  Some  of  the  many  decision 
problems  that  might  arise  in  planning  and  designing  service  contracts  are: 

i)  Design  of  the  contract. 

ii)  Capacity  planning  for  repair  and  service  facilities. 

iii)  Inventory  policies  for  replacement  parts  and  repair  resources. 

iv)  Optimal  cost-reliability  tradeoff. 

v)  Optimal  product  testing  and  quality  control  effort. 

vi)  Design  of  data  collection  and  information  systems. 

vii)  Special  problems. 

In  making  these  decisions  the  manufacturer  should  employ  all  the  information 
he  can  obtain  about  the  failure  characteristics  of  the  product.  If  it 
possible  to  model  and  predict  part  failure,  the  relationship  between 
product  cost  and  reliability,  expected  costs  of  repair  etc.,  the  decisions 
involved  can  be  made  on  a  systematic  basis. 

DECISION  PROBLEMS  IN  PLANNING  AND  DESIGN: 


All  the  decisions  described  below  can  be  seen  to  be  highly  interrelated. 
However  the  issue  of  interactions  is  deferred  for  the  present  for  the  sake 
of  expositional  clarity. 

1)  Design  of  the  Service  Contract  : 

This  relates  to  the  design  of  the  terms  of  service  contract  offer. The 
main  variables  involved  are  : 

i)  Price 

ii)  Time  horizon 

iii)  Extent  of  services  to  be  provided. 
Price  :      This  is  naturally  related  to  the  other  terms  of  the  contract. 
Once  those  are  fixed,  the  price  must  be  set  not  only  on  the  basis  of  profit 
margin  considerations  and  cost  of  providing  the  services  contracted  for, 
but  also  on  the  marketing  and  promotional  value  of  the  contract  and  other 
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competitive  factors.  Thus  the  setting  of  the  contract  is  a  highly  complex 
decision  involving  many  subjective  factors.  However  important  inputs  to 
the  decision  are  the  cost  of  providing  services,  and  the  demand  curve  for 
the  service  contract. 

Time  Horizon  :       Clearly,  the  contract  cannot  cover  all  servicing 
to  infinity  and  in  general  the  longer  the  time  of  coverage,  the  greater 
the  marginal  cost  to  the  manufacturer  and  the  greater  the  price  of  the 
contract.  It  would  probably  be  necessary  to  cut  off  coverage  before 
failure  rates  of  the  product  start  to  increase  drastically,  and  the  time 
of  coverage  is  thus  also  informative  to  the  consumer.  There  are  many  other 
factors  that  intervene;  for  example  even  if  the  product  is  discontinued, 
service  and  spare  parts  have  to  be  made  available.  From  the  point  of 
view  of  consumer  perception  and  behavior,  it  may  only  be  necessary  to 
ensure  coverage  for  the  first  couple  of  years  after  purchase  to  provide 
the  buyer  with  enough  security.  This  is  especially  true  of  durables  which 
are  highly  vogue  dependent  and  have  high  rates  of  obsolescence.  The  key 
inputs  to  this  decision  would  be  failure  models,  the  costs  of  service  for 
different  time  horizons,  models  of  consumer  behavior  in  terms  of  original 
buyer  usage  and  total  usage,  etc. 

Extent  of  Services  :     Coverage  should  presumably  be  sufficiently  broad 
to  remove  from  the  buyer  any  sense  of  insecurity  or  uncertainty;  i.e.  the 
buyer  should  be  protected  from  the  onus  of  bearing  large  costs  due  to  random 
failures.  However,  the  time  of  coverage  might  be  set  differently  for  different 
parts  and  subsystems,  and  some  items  could  conceivably  be  excluded  from 
coverage.  Again  failure  rate  models  are  important  in  making  these  decisions. 

2)  Facilities  Planning  : 

Facilities  supporting  a  service  contract  offering  would  include  repair  shops, 
spare  part  and  repair  resource  inventories,  product  collection  and  return 
systems,  manpower  organization,  data  collection  and  information  systems  etc. 
The  capacities  required  are  determined  based  on  the  anticipated  volume  of 
requests  for  repair  and  service  under  the  contract.  These  in  turn  depend 
on  the  number  of  contracts  in  force  and  rates  of  part  and  subsystem  failure. 
We  thus  need  models  of  demand  for  the  contract  as  well  as  failure  models. 


3)  Inventory  Policies  : 

In  determining  inventory  policies,  the  demand  for  replacement  parts  and 
material  requirements  can  be  predicted  from  knowledge  of  failure  charac- 
teristics. Given  the  population  of  units  in  service  the  demand  for 
service  can  be  predicted  based  on  statistics  from  failure  models  such 
as  renewal  or  restoration  processes. 

4)  Production  Cost  -  Reliability  Tradeoff  : 

The  less  reliable  a  product,  the  greater  the  costs  to  the  manufacturer  of 
providing  service  and  repair.  If  the  price  and  time  horizon  of  the  service 
contract  are  determined  by  considerations  other  than  these  costs,  then 
it  is  worthwhile  to  build  greater  reliability  into  the  product.  To  decide 
on  the  appropriate  level  of  effort  in  this  area,  a  model  is  required  that 
relates  reliability  to  design  and  production  costs. 

5)  Product  Testing  and  Quality  Control: 

Testing  provides  improved  knowledge  about  the  reliability  of  the  product. 
In  general  better  knowledge  and  reduced  uncertainty  reduce  costs  by  allowing 
planning  to  be  done  more  precisely  with  less  costs  of  safety  margins  and 
hedging.  To  cite  a  simple  example  —  better  predictions  of  demand  for 
replacement  parts  reduce:safety  stock  requirements  and  hence  the  associated 
holding  costs.  Similarly,  improved  quality  control  procedures  would  lead 
to  a  reduction  in  service  and  repair  requirements.  In  formulating  decision 
problems  in  these  areas,  we  need  models  that  enable  us  to  place  a  value 
on  information  about  failure  rates. 

6)  Data  Collection  and  Information  Systems: 

Incoming  consumer  requests  for  repair  and  service  are  informative  in  them- 
selves and  ought  to  be  tracked.  Here  it  is  necessary  to  determine  what 
data  is  significant  and  how  it  is  to  be  used  in  improving  decision  making. 

7)  Special  Problems  : 

Unforeseen  situations  with  large  associated  losses  may  suddenly  occur  as 
in  the  case  of  the  detection  of  a  hazardous  defect  in  the  product.  These 
cannot  be  explicitly  planned  for  except  from  the  viewpoint  of  maintaining 
a  tracking  and  detection  facility  included  in  the  data  collection  system. 


THE  MODELLING  APPROACH  : 

The  preceding  discussion  of  decision  problems  has  revealed  a  need  for 
models  in  the  context  of  which  decision  problems  can  be  framed. Here  we 
we  discuss  an  approach  to  model  building  in  terms  of  the  basic  models 
that  appear  to  be  required  . 

1)  Failure  Rate  Models  : 

We  require  descriptive  models  of  failure  of  parts  in  the  form  of  probability 
distributions  of  the  time  to  failure  (TTF)  of  each  product  element.  There 
may  be  dependencies  between  failure  of  different  parts  and  joint  distri- 
butions may  be  appropriate.  Depending  on  the  terms  of  the  service  contract, 
we  may  need  models  of  subsystem  failure, "catastrophic"  or  final  failure, 
etc.  In  addition  by  using  Bayesian  methods,  we  can  incorporate  the  concept 
of  updating  astate  of  knowledge  as  information  becomes  available. 

2)  Sequences  of  failures  : 

The  occurrence  of  failures  of  a  part  followed  by  repair  or  replacement 
gives  rise  to  a  sequence  through  time  which  can  be  modelled  as  a  renewal 
or  a  restoration  process.  Again  a  Bayesian  treatment  would  allow  us  to 
update  our  knowledge  of  the  process  as  we  observe  it.  Statistics  from 
these  models  can  be  used  to  predict  aggregate  demand  for  service. 

3)  Sequences  of  costs  : 

Each  failure  and  service  request  gives  rise  to  an  associated  cost  of  repair 
or  replacement  [r/r)   The  present  value  of  such  a  sequence  of  future  costs  is 
an  important  input  to  various  decisions.  The  expected  value  of  these  cousts 
discounted  to  the  present  is  a  very  basic  model  in  the  present  context. 

4)  Production  Cost/Reliability  Model 

One  way  to  approach  the  problem  to  relating  costs  of  production  and  design 
to  product  reliability  might  be  through  the  prior  distributions  for  the 
parameters  of  the  time  to  failure  (TTF)  model.  We  could  subjectively 
estimate  the  prior  conditional  on  the  production  costs;  or  we  could  assert 
that  increased  production  and  design  efforts  lead  to  some  change  in  the 
prior  parameters.  Alternatively  we  might  argue  that  over  a  small  range, 
the  time  to  failure  changes  linearly  with  incremental  change  in  design  and 
material  costs,  and  use  Bayesian  regression  models.  However  this  type  of 
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model  would  greatly  complicate  the  analysis,  and  furthermore  in  the 
interests  of  tractability  we  would  have  to  restrict  considerations  to 
Normal  probability  distributions. 

5)  Value  of  Sampling  and  Information  : 

Bayesian  models  enable  us  to  treat  questions  such  as  the  value  of  information, 
gain  from  sampling,  and  optimal  sampling.  These  models  are  the  basis  for 
tradeoffs  involving  costs  of  information  gathering  as  in  the  cases  of 
testing,  quality  control  and  data  collection  that  were  described  earlier. 
The  value  of  information  is  determined  by  comparing  expected  cost  without 
the  information  against  expected  cost  under  the  updated  state  of  knowledge 
after  obtaining  the  information. 

6)  Demand  For  Service  Contracts  : 

The  demand  for  service  contracts  determines  the  overall  capacity  levels 
required  in  planning  for  service  and  repair  facilities.  This  demand  is 
a  function  of  the  price  of  the  contract,  and  consumer  perceptions  of 
reliability  and  repair  costs.  It  is  clearly  bounded  from  above  by  the 
demand  for  the  product  itself.  Regression  models  might  be  useful  in 
prediction  of  this  demand  based  on  variables  such  as  price,  advertising 
and  promotional  efforts,  rate  of  obsolescence  etc. 


SPECIFIC  MODELS 

In  this  section  we  describe  some  modelling  efforts  in  greater  detail. 
The  basic  models  relate  to  part  failure  and  sequences  of  failures  with 
their  associated  costs. 

Time  to  Failure  (TTF)  Models  : 

The  basic  descriptive  model  of  part  failure  consists  of  a  probability 
distribution  on  the  time  to  failure  (TTF).  The  literature  in  this  area 
is  extensive  and  a  wide  variety  of  distributions  has  been  used  in  modelling 
failure  characteristics. (  Barlow  and  Proschan[5],  Zelen  [6]).  The  Expo- 
nential, Erlang,  Weibull  and  Normal  distributions  are  commonly  employed 
failure  models. 

For  our  present  purposes  we  are  interested  in  Bayesian  models  of  failure. 
This  in  effect  involves  considering  the  parameters  of  the  TTF  distribution 
as  not  known  with  certainty.  Prior  distributions  are  assigned  to  the  un- 
known parameters  and  are  updated  as  new  information  becomes  available. 
Models  such  as  the  exponential  and  Normal  distributions  are  well  suited 
to  Bayesian  modification.  The  Weibull  distribution  however, does  not 
possess  sufficient  statistics  of  fixed  dimensionality,  and  hence  is  not 
amenable  to  the  usual  form  of  natural  conjugate  analysis  although  Bayesian 
treatments  are  available  in  the  literature  (Sol and  [7],  Bassin  [8]). 

The  type  of  model  chosen  in  a  particular  case  will  naturally  depend  on 
the  failure  characteristics  of  the  product  in  question.  However  we  are 
here  especially  concerned  with  the  early  stages  of  product  lifetime,  and 
this  emphasis  warrants  some  discussion  in  terms  of  its  effects  on  model 
selection.  Most  products  are  subject  to  "infant  mortality";  i.e.  there 
is  a  segment  of  the  population  of  all  units  produced  that  fails  very 
early,  possibly  due  to  some  undetected  defect.  In  the  case  of  complex 
products  with  many  parts  there  may  be  many  distinct  part  failures  in  a 
short  time,  or  there  may  be  persistent  failure  of  a  particular  part  or 
subsystem.  In  other  words  the  population  of  all  products  contains  some 
.raction  of  defectives  or  "lemons". 


The  effect  of  such  a  subpopulation  of  defectives  on  failures  observed  is 
shown  in  Figure  2  which  shows  a  record  of  cumulative  service  requests(in 
terms  of  %of  total  sales)  versus  time  measured  from  date  of  sale. The  num- 
bers are  hypothetical  but  the  general  characteristics  of  the  curve  are 
fairly  typical. 
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Figure  2 

The  initial  transient  portion  of  the  curve  suggests  the  existence  of 
defectives  with  a  relatively  high  failure  rate.  As  units  are  repaired 
and  returned  to  the  population,  the  proportion  of  defectives  and  the 
overall  failure  rate  both  decrease.  The  asymptotic  straight  line  portion 
would  thus  represent  a  constant  average  failure  rate  which  may  be  inter- 
preted as  being  due  to  the"usual"  incidence  of  random  failures. 

The  Weibull  process  is  often  used  to  model  such  behavior  exhibiting 
decreasing  failure  rates.  However,  it  is  not  an  analytically  tractable 
model,  especially  in  the  Bayesian  context,  and  estimation  of  its  para- 
meters is  a  cumbersome  process  (  Hausman  and  Kamins  [16]). In  looking 
for  alternative  models,  it  was  hypothesised  that  a  failure  distribution 
consisting  of  a  mixture  of  exponential  densities  would  exhibit  the 
desired  failure  rate  characteristics  and  this  can  easily  be  verified. 
Furthermore  we  can  interpret  this  model  intuitively  as  representing  a 
mixture  of  two  kinds  of  units  --  a  proportion  'p'  of  defectives  with  a 
high  failure  rate  (A,)  and  a  proportion  '1-p'  of  "normal"  units  with  a 
lower  failure  rate  (A2)  ascribed  to  random  failures. 
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Sequences  of  Failures  and  Costs: 

The  sequence  of  failures  followed  by  repair  or  replacement  so  that  the 
unit  is  "as-good-as-new"  constitutes  a  renewal  process.  This  type  of 
process  has  been  extensively  modelled  and  analysed  in  the  literature 
(Cox  [9],  Zelen  [6]).  We  may  also  wish  to  consider  restoration  process 
models  where  servicing  does  not  leave  the  unit  as  good  as  new  but  rather 
restores  it  to  its  state  at  the  time  of  failure  so  that  it  is  "as-bad- 
as-old"  (Ascher  [10]).  Suppose  that  failure  of  a  unit  can  be  caused  by 
one  of  a  number  of  parts.  Then  failure  of  one  part  followed  by  its  rep- 
lacement does  not  render  the  unit  as  good  as  new  and  hence  a  restoration 
model  might  be  appropriate.  Both  types  of  models  have  been  studied  in 
a  Bayesian  context  (Bassin  [11]).  In  the  following  however,  we  shall 
chiefly  be  concerned  with  the  renewal  model. 

Consider  the  stochastic  failure  of  an  element  of  a  system.  Suppose  that 
each  failure  is  followed  by  (instantaneous)  repair  or  replacement  of  the 
element.  Then  the  resulting  sequence  can  be  modelled  as  a  renewal  process 
with  associated  costs  of  renewal.  Many  of  the  results  that  follow  are 
essentially  available  in  Cox[9]  but  are  repeated  for  continuity;  the 
notation  we  adopt  paraphrases  that  used  by  Cox  : 

t;.  :  time  to  failure  of  the  i'th  component. 
ft*):  p.d.f.  (of  time  to  failure) 
I£t|e):  p.d.f.  conditional  on  parameter  set 
<r'(9le>'):  Prior  distribution  on  &    in  the  Bayesian  sense,  with 

prior  parameter  set  _£* 
D(t[^1):  p.d.f.  of  £  unconditional  on  the  uncertain  parameter 

set 
F(t)  :  c.d.f.  (of  time  to  failure) 

G^)  =  jf  Ki)<K  -  PLt»t3-i-Fft) 

2C±)  •    (Instantaneous)  Failure  rate  or  Hazard  rate  function. 

<fo)  =  d{$(');  $\     •   The  Laplace  transform  of  a  function  $0) 
?<*>-  Etc"**]  -  ±{  K-);s] 

s 


S-,  ■   Total  time  to  the  r'th  renewal 
S,  -  t,  +  t4+  ..  .  +  tn 

fe^t)  :.p.d.f.  of  ?, 

KY(t)=  c.d.f.  of  S, 

fjfc  :  The  number  of  renewals  in  (0,t) 

H(t)  •  The  renewal  function:-  the  expected  no. of  renewals  in(0,t) 

-UCt)  =  i.  H(t)  :  the  renewal  density 

of 

L(»  =  E(Nt2  )  1  VH(Nt)s  L(t)-  LHCt)]"*" 

With  regard  to  the  sequence  of  costs  associated  with  renewals,  we  are 
interested  in  the  expectation  of  future  costs  incurred,  discounted  to 
the  present.  There  are  several  different  quantities  that  are  of  interest. 
We  may  want  to  know  the  cost  of  an  infinite  sequence  of  renewals  or  of 
a  finite  number  of  renewals.  In  the  service  contract  context  in  particu- 
lar, we  might  want  to  know  the  discounted  costs  associated  with  a  sequence 
of  renewals  that  is  truncated  at  a  finite  time,  without  restriction  on 
the  number  of  renewals.  Let  us  denote  : 

R(t)  :  the  cost  associated  with  a  renewal  at  time  t.  We  will 
assume  for  the  present  that  this  is  constant,  R£t)=  R 

^  :  discounting  factor  for  continuous  discounting  of  future 
costs. 

"V-Y  :  discounted  cost  of  the  r'th  renewal; 
V-<  :  total  discounted  cost  of  the  first  r  renewals 

V   :  TDC  of  an  infinite  sequence  of  renewals 

-    Z^ 

\y(T)  :  TDC  of  a  sequence  of  renewals  truncated  at  time  T (but 
with  no  limit  on  the  number  of  renewals). 

As  usual  the  tilde  is  used  to  identify  random  variables  (r.v.'s).  A  bar  will 
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be  used  to  denote  expectation  of  a  r.v.:  -vY  represents  the  expected 
discounted  cost  of  the  r'th  renewal,  V   is  the  EDC  of  a  sequence 
truncated  at  time  T  and  so  on.  We  shall  assume  that  the  TTF's  of 
successive  elements  are  mutually  independent,  identically  distributed 
(m.i.i.d.)  r.v.'s,  conditional  on  the  parameters  of  the  distribution 
being  known.  We  first  derive  general  formulae  for  the  quantities  of 
interest,  and  then  apply  them  to  particular  failure  distributions. 

CASE  I  :  Parameters  of  the  TTF  density  known  with  certainty 

For  a  single  failure  followed  by  renewal  : 

■v-    ,     R  i*tf)  0) 

r  /  X   ' 

Since  the  t^s  are  m.i.i.d.  r.v.'s  the  transform  of  the  sum  of  the  Xc  s 
is  equal  to  the  product  of  the  transforms,  and  hence  : 

^  =  fcjfc*  Cf )  =  R.TT  i*(?) 

For  a  sequence  of  n  renewals,  the  total  expected  discounted  cost  : 
and  hence  for  an  infinite  sequence  of  renewals, 

i-P<»  _ft 

We  note  that  since  for  the  usual   values  of  f  t  e     $(*" )  <  ££fc)  ,  we  have 
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To  find  the  variance  of  this  quantity  first  note  that  : 

For  j>(.  ,  we  can  write  Sj  *  S;  + 1^+  +  ..  .  + 1  •  /  and  hence  the  second 
term  in  the  equation  above  can  be  replaced  by  X  2.  '2-     <~  £    w"' 

Whence  applying  the  expectation  operator 

which  upon  simplification  gives  : 


E(v-).   RM*(in 


\  -pr? )  J 

and  the  variance  can  be  found  in  the  usual  way. 

If  we  now  consider  the  sequence  of  renewals  truncated  at  time  T,  then 

■1    Z.\       O 

which  can  also  be  written  as 

VfT)  ,  fcjVftJk(t-)  =it  On 

0 

CO  , 

where   J^(t\   "S"  U  /(-)  =  °L  m^t  )   "»  the  renewal  density. Clearly 

we  could  easily  replace  ^  with  some  time  dependent  function  R-Ct). 
Again  in  determining  the  variance  of  this  quantity,  since  by  truncation 
we  are  saying  that  the  cost  of  a  renewal  for  t>T  is  zero  we  can  write: 

"  Ul    0  I-  I  1  -  «-+  '   0  g 

which  on  simplification  yields 

fltM'l.      ft*[i-Kn1,  fY2<?frA(t)jr  rt) 

L  J  "  [»  -l(f)l    i 

where     : 


Considering  now  the  number  of  renewals  in   (0,t)  :    (Cox[9]) 

We  know  that   Nt<.T  implies  S^>t    .  Therefore  ^[Nt<v]=  P  [S-r>  t]--  I  -  XyLt) 

and  hence  we  can  write  P[wt=>]r  Kytt)-  k      (t)  -We  can  thus  determine 

the  renewal    function   :  ^ 

Hit),    ECNt)-.   £^P[Nt  =  ^]    =  £^I>vtO-Kv+^)] 

WCt)--     £    K    (t  )  (»o) 


Taking  transforms  on  both  sides   : 
* 


rn) 


H^0  s  LIAVO 

and  we  know  that  -fe*Y  (O- [  £*(*•)]   hence 

sli  -  iWD 
HCt)=  f>no] 

In  finding  the  variance  of  the  number  of  renewals,  first  find 
L(t)  =  E(V)  =  Z^V[^vCt)-WY+,:'t)] 

■V  s  0 

Taking  transforms  ,  simplifying  and  noting  that    ZL   "*  P   ~   P  ^  +  ^  - 

v:o  It  -  P  )  3 

we  get 

■)- 


and  we  can  thus  find  the  variance  as 


Van 


t5t).  uo-  |>>r     ;       uo.i-'^'l 
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The  distribution  of  the  number  of  renewals  may  also  be  obtained  by 
studied  the  generating  function  of  Nt  (  Cox[9]).  The  results  are 
cumbersome,  and  we  will  not  generally  obtain  this  distribution.  How- 
ever, we  briefly  indicate  the  procedure  involved.  Let  us  denote 


v:  \ 

Taking  transforms  on  both  sides  : 

-   i  -  j»(o 

And  we  can  find  PlN^J  as  the  coefficient  of  £   in   H£t,  2  ) 

CASE  II  :  Parameters  of  the  TTF  distribution  not  known  with  certainty 

One  approach  is  to  assign  a  prior  distribution*  (©.  I  fi' )  to  6_  the  set 
of  parameters  not  known  with  certainty,  and  then  to  derive  the  uncondi- 
tional distribution  of  time  to  failure  for  a  single  failure. 

® 


so  that   — 


v(  =  RD*«  )  (it) 

For  a  sequence  of  failures,  the  unconditional  distributions  of  the 
J.otal  times  to  failure  for  successive  failures  are  not  independent,  since 
they  share  the  same  prior  distribution  on  the  uncertain  parameter  set. 
Hence  the  simple  results  obtained  in  the  previous  case  do  not  apply  when 
we  use  this  direct  approach: 
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VT  ^  R  £  D?(f)  08) 

i  =  i 

However  a  simple  alternative  approach  is  just  to  take  the  prior  expectations 
of  the  expressions  obtained  conditional  on  the  parameter  being  known  (Case 
I).  Thus  for  example: 

({  }«e-**,(t|g)it}  rCfi|i»J)Je 


(9    o 


And  similarly  : 


1  V  *vf,n  I 


vn)  s   E  [v(t;(q)  ] 


f22) 


When  considering  variances  of  quantities,  we  have  to  use  the  appropriate 
formula  for  iterated  variances. 

Applications  to  particular  TTF  models 
(A)  The  EXPONENTIAL  TTF  distribution: 

At-  \  "At 
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For  the  exponential   density: 

■st-    -      ~xfc  ■■ 


0 

and  hence  direct  substitution  in  equations  (1),(3),(4)  and  (5)  gives: 


-v,  =   RX 


v,.   *Mi 


■v; 


■  ■&>)■ 


From  equation  (7)  we  also  have 

and  hence  the  variance  of  the  total  discounted  cost  of  an  infinite 
sequence  of  renewals 

"■'"        ~        Tr 

To  compute  the  quantity  V   we  first  find  the  renewal  density  : 

»   ,  p.  ,.L>i.|  -At 

^(t)  -  £*-,(*;    =  ^  (>-D i.  a 

And  by  equation  (8)  : 

o 

We  can  find  the  variance  of  this  cost  by  using  equation  (9)  which  gives 


F_(v^2]  =  HA(^-e"2?T)  ['  +*<«)  1 

2-r         y  Li  -i/r )  J     c 


*9) 


where  :    I(r\  B  A  ( i  —  sT'*"1  f)"1"  \ 

The  variance  can  now  be  found  in  the  usual  way. 

In  determining  the  renewal  function,  we  note  that  we  have  already  obtained 
the  renewal  density  k(t  )  ■=  X   ,  and  hence 

To  find  the  variance  of  Nt  ,  equation  (12)  gives 
Therefore 


Suppose  now  that  the  parameter  ^  is  not  known  with  certainty,  and  we 
assign  to  it  a  Gamma-1  natural  conjugate  prior  density  .(Raiffa  & 
Schlaifer  [12]): 

_Xt'  ,i  , 

5\  -fyiC^W.t')  OC  *.  *  ;      *,V,t'>0      (32.) 

The  unconditional  density  of  the  total  time  to  the  r'th  renewal  is 
given  by  : 

DvOKt';v)=    JfwCth?0X(*Kt,W>  C33) 

o 

and  has  been  shown  in  [12](7.4.2)  to  be  of  the"inverted  Beta-2"  form: 


^CtKt';v)=  ft.aCt|Tfv#t')-  ■ 


_j t'  '  t' 

BO,  V)  0+fJy+Y 


We  could  presumably  obtain  directly  quantities  such  as  : 

o 

However  this  "direct  approach"  is  explored  in  Appendix  A  and  is  found 
to  be  quite  cumbersome.  Noting  that  the  parameter  X  appears  in  a  simple 
manner  in  equations  (25)  -  (31)  we  can  simply  write  : 


v  ■  M£)  ■ 


s 

where  >'=  j'    is  the  prior  expectation  of  X  .  Other  results  follow 

T1 

in  a  similar  manner.   For  example  : 


*(££) +  V->(V) 


^«*(>/)   =        ^N'    f  z+-  -ft')  r3?' 
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(B)  The  MIXED  EXPONENTIAL  TTF  distribution: 

This  distribution  was  suggested  as  a  model  adequate  to  describe  failures 
exhibiting  a  decreasing  failure  rate.  As  stated  earlier,  the  intuitive 
notion  motivating  this  model  conceives  of  the  population  of  all  parts  as 
being  composed  of  two  sub-populations  :  a  fraction  p  subject  to  both 
defective  and  random  failure,  and  a  fraction  (1-p)  subject  only  to 
random  failure.  If  it  is  assumed  that  both  sub-populations  exhibit  expo- 
nential failure  distributions  with  rates  A ,  and  >\x   respectively,  we 
may  write  : 

Jft)  =  F4,(t)  +(>-p)j,(t)  =  pA,c"A,t  *0-P)Axe"  x 

We  may  generally  assume  that  *,>*2.  and   p<  f »  -  p.) 

We  will  first  examine  the  failure  rate  characteristics  of  this  distribution 

At  t-.o  we  see  that  ££o)r  p  x,+  6-  p  )>  j_  too) 

Furthermore,  Z(t)  is  monotone  decreasing  in  t  since  considering  the  sign  of 

the  derivative  of  Z(t)  : 

-f\+Ax)t 

^-2Ct)   =        -  CA.-Q  f  Q-P)e <       0  (41) 

[  Pt    '     +  (>-p)c   a    J 

Rewriting  this  as 

2L2Ct  )    =         -(\-^'?(|-f) 


at  iP0-p)<  p=-e-K-Ot  +(,.piV(->"-x'>t 


We  can  see  that  for  any  X,,  Xj. 
and  we  can  also  see  by  writing 


0  (4») 
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Ut) 


PA'*         +  Q-f)*a. 


that  for   >,  >  >2_ 

t  -*<» 
We  define  here  for  future  convenience  the  quantity 

z. 


D 


(«3) 


f4<0 


4<n 


t  =  o 


which  might  be  interpreted  as  a  sort  of  measure  of  the  heterogeneity  of 
the  population  being  modelled. 

We  have  thus  shown  in  equations  (40)  through  (43)  that  the  mixed  exponen- 
tial failure  density  exhibits  a  decreasing  failure  rate  : 


Figure  3 


Let  us  now  consider  the  behavior  of  the  renewal  function  for  this  density: 


res)  =  _L±L  *  L'-P)>- 


Using  equation  (11 ) 


whence  after  some  algebra  we  obtain 

W\S)    =       >■>,.  +     sl>V  (Up)Xj 


'M) 


(Ml) 


-21- 


Defining    -A,  =  P^+(l"P)Xi  (A8) 

7^  .  (i-p)A,  +  PX,  =  (>,  *M~-A-, 
We  can  rewrite 

To  take  the  inverse  transform,  we  first  note  that 


rso) 


and  hence  upon  simplification  we  can  write  the  renewal  function  as 
Investigating  the  coefficient  of  the  exponential  decay  term  : 


(**) 


The  renewal  density  is  given  by  the  slope  of  the  renewal  function: 
and  we  can  see  that 

-k(o)  --  V\_v  =  2^) 

(54) 

The  shape  of  the  renewal  function  thus  looks  yery   much  like  the  curve 
in  Figure  2  in  its  general  behavior,  and  is  shown  in  Figure  4.  Whether 
Jiis  model  will  in  fact  fit  actual  data  remains  to  be  verified. 
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H(t) 


Figure  4 


In  computing  the  variance  of  the  number  of  renewals,  substituting  the 
appropriate  quantities  in  equation  (12)  we  get  after  some  algebra  that 

L^  "    (7~>  {  A  ■  +  l^M^i)  ] 


(SS) 


To  write  the  inverse  transform  let 


and  note  that 


fsM 


•Axt 


a2tt)  = 


t«£ 


Axt 


■A, 


t") 


^  ■  i-  l*tol-  £-  V1*^*^0' 
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whence  we  can  write 

L(t)  =   Z_  X.  «-'Jt)  (58) 

I  =  \ 

Examining  now  the  cost  models  corresponding  to  this  failure  density,  we 
can  write  from  equations  (5)  and  (8)  that 


We  will  not  write  out  the  variances  of  these  quantities  in  all  detail. 
We  remark  however  that  in  equation  (7)  the  multiplying  factor  can  be 
written  as 


|  LLC 


and  in  equation  (9)  we  can  write  for  the  integral 

OXOt  i    f         a  r,    -  Cf-tX^T 


\  -  ret) 


Parameters  not  known  with  certainty  : 

Estimation  of  the  parameters  for  the  mixed  exponential  density  is  a  non- 
trivial  task.  While  this  distribution  possesses  desirable  characteristics 
and  gives  a  more  intuitive  interpretation  of  the  failure  process  than 
say  the  Wei  bull  distribution,  it  does  not  buy  us  tractability  in  a  Bayesian 
approach.  The  distribution  does  not  possess  sufficient  statistics  of 
fixed  dimensionality  and  hence  cannot  be  assigned  a  natural  conjugate 
prdor  in  the  usual  way.  However,  Behboodian  [17]  has  investigated  the 
Bayesian  treatment  of  the  mixing  parameter  p  and  shows  that  the  use 
of  a  Beta  prior  on  p  leads  to  a  mixed  Beta  posterior  density.  Thus  the 
ramily  of  mixed  Beta  densities  is  conjugate  to  the  likelihood  function 
for  the  mixing  parameter  p  of  a  mixture  of  distributions.  The  number  of 
terms  in  the  posterior  mixture  depends  on  the  sample  size. 
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The  estimation  problem  in  a  classical  setting  has  been  studied  by  Boes 
in  [18]  and  [19].  Rider  [20]  has  examined  the  application  of  the  method 
of  moments  and  Mendenhall  and  Hader  [21]  have  studied  the  maximum  like- 
lihood approach  to  the  mixed  exponential  case. 

In  general  we  may  obtain  information  about  the  process  in  different  ways. 
For  example  the  data  could  be  from  : 

i)  A  random  sample  that  is  subjected  to  a  life  test  for  a  fixed  period 
at  the  end  of  which  there  may  or  may  not  be  survivors, 
ii)  Historical  data  on  the  aggregate  rate  of  demand  for  service  requests; 

i.e.  cumulative  failures  as  in  figure  2. 
iii)  Presently  incoming  data  on  service  requests,  where  it  may  be  possible 
to  obtain  some  additional  information  on  the  failure. 

The  treatment  of  the  problem  varies  depending  on  the  particular  situation, 
and  the  assumptions  made  about  the  sampling  process.  For  example  if  we 
assume  that  a  failure  can  always  be  identified  as  having  occurred  from 
one  or  the  other  of  the  subpopulations,  the  problem  is  considerably 
simplified  and  becomes  a  routine  case  of  prior  to  posterior  analysis. This 
assumption  may  not  however  be  very  realistic.  In  particular  we  note  that 
for  the  interpretation  suggested  here,  it  is  possible  for  a  defective 
element  to  fail  in  a  random  mode  and  be  indistiguishable  from  a  random 
failure  although  the  reverse  might  not  be  the  case. 

The  practical  importance  of  estimation  of  the  parameters  is  directly 
related  to  the  decisions  that  are  to  be  made.  It  may  be  that  A  2.  the 
failure  rate  for  random  failures  is  quite  well  known  from  experience 
with  similar  products.  In  any  case  it  is  probably  not  a  controllable 
variable.  The  parameters  A,  and  p  on  the  other  hand  might  be  considered 
control  variables  in  the  following  sense  :  If  more  effort  is  put  into 
testing  screening  and  quality  control  one  could  hope  to  reduce  the 
fraction  of  defective  items  in  the  population.  Similarly  the  control  of 

>  1  might  be  achieved  as  a  result  of  deliberate  manufacturing  policy 
interms  of  reliability  and  durability  designed  into  the  product.  Further- 
more testing  can  improve  knowledge  about  these  parameters  and  this  will 
affect  the  decisions  that  are  made.  Thus  we  see  that  these  parameters  have 
a  natural  interpretation  in  a  decision  making  context. 
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MODIFICATIONS  &  EXTENSIONS 


We  have  outlined  the  basic  approach  to  modelling  sequences  of  failures 
and  costs  with  a  view  to  formulating  decisions  about  planning  and 
designing  service  contract  offerings.  It  is  relatively  simple  to  make 
certain  modifications  to  these  models  to  fit  particular  situations. 
Some  of  the  many  possibilities  are  briefly  sketched  out  below  : 

a)  Service  requests  and  costs  for  a  group  of  identical  products  : 
Assuming  that  these  fail  independently,  the  expected  number  of  requests 
for  service  is  given  by  the  sum  of  the  expectations.  Similar  results 
hold  for  the  variance  of  the  number  of  requests  ,  expected  costs  and 
variance  of  the  costs. 

b)  Different  dates  of  sale  : 

If  the  sale  dates  are  known  for  two  batches  say,  the  expected  number  of 
requests  can  be  found  simply  by  adding  graphs  such  as  in  figure  4,  suitably 
translated  along  the  time  axis  to  account  for  the  date  of  sale.  If  the 
sale  date  is  a  random  variable,  then  we  can  use  the  artifice  of  regarding 
sale  of  the  product  as  the  first  "renewal".  The  process  can  now  be 
modelled  as  a  "modified"  renewal  process  (Cox  [9]  ).  The  date  of  sale  is 
in  a  sense  a  control  variable,  since  by  staggering  sales  in  time,  it  is 
possible  to  smooth  out  the  demand  for  service  requests. 

c)  Different  failure  densities  for  original  and  replacement  parts  : 
Suppose  we  can  assume  that  a  repaired  product  has  a  different  failure 
density  from  the  original  --  for  example  we  might  assume  that  there 
are  no  defectives  amongst  replacement  parts.  Again,  a  modified  renewal 
process  would  be  an  appropriate  model. 

d)  Different  cost  functions  : 

We  can  see  from  equations  (8)  and  (28)  that  it  is  fairly  straightforward 
to  consider  other  than  constant  costs.  In  particular  we  may  want  to  use 
polynomials  or  exponential  functions  in  t  to  reflect  costs  which  change 
systematically  over  time.  We  may  also  wish  to  consider  stochastic  costs 
ind  utilities. 

Apart  from  these  modifications,  it  may  be  of  interest  to  determine  other 
information  from  the  same  basic  model.  For  example,  the  average  age  of 
parts  in  service  can  be  determined  by  finding  the  mean  "backward  re- 
currence time  ". 
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APPENDIX  A  : 

In  this  appendix  we  explore  in  some  detail  the  approach  suggested  in 
equation  (35).  That  is  to  say  we  attempt  to  obtain  some  expressions 
for  future  discounted  costs  for  the  case  of  exponential  TTF  distri- 
butions, by  working  directly  with  the  unconditional  sampling  density 
obtained  by  assigning  Gamma-1  priors.  The  failure  density  unconditional 
on  the  parameter  A  is  of  the  "inverted  Beta  2"  form.  Recalling 
equation  (35)  : 


Let  us  denote  the  integral:   ^ 

/  -S  x   "»-\     i  / 

O       (X  -V  o   ')~',  T 

so  that 

For  integral  values  of  r  and  r'  Erdelyi  et.al.[14]  give  that 

Now  r  is  integer  but  in  general  r'  is  not.  In  this  case  we  obtain  from 
Gradshteyn  and  Ryzhik  [13]  that 

!(■*,■*•  o  =  lc^'i;l  c^ITo  W,.*v->'  a.'^i>)        ;a.v 

where  W  0)is  the  Whittaker  Function  which  is  written  as 


where 


The  confluent  hypergeometric  function  ,F,  denotes  the  series 

,F  io/  V   ?  ^  -   I  4-   *  *  +  ±j£±0  *%  ^(c,tl)^+;)Z  '  +  ..    U,7) 

Substituting  these  expressions  in  A  A  we  have  : 


^f'V.c*,.-*';  *,)>*]     i*;.; 
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For  the  case  r=l  we  can  directly  evalute  the  integral 

v;  .-   Uest  -J ill   ,  AC 

i     &ov )  (t  4  t')i+y 


T(- 


which  can  be  found  from  the  handbooks  to  be 

1  (:/^"  -  iif/^'Etf-ft1;    ca.»o3 


where  the  exponential  integral  is  denoted  by 


-Ed-It')   -  E,(fi'>  ■   J  t«"V  ^         ^U) 

■ft' 

Clearly  evaluation  of  V  say  by  direct  summation  of  costs  is  going  to 
be  quite  involved  due  to  the  extreme  unwieldiness  of  these  expressions. 
Of  course  we  have  already  determined  V  in  equation  (36).  We  shall 
here  demonstrate  that  we  can  reproduce  these  simple  results  by  directly 
working  with  the  inverted-Beta-2  density.  Let  us  write 

Co    <x>  I 

Reversing  the  order  of  summation  and  integration  and  simplifying 

where  we  have  again  set  a  =  t'  for  notational  convenience. 
Now  this  interchange  is  feasible  if  and  only  if  the  series  in  brackets  is 
uniformly  convergent  in  the  the  sense  that  the  sequence  of  functions  formed 
by  the  partial  sums  of  the  series  converges  uniformly  (Rudin  [15]}.  In 
fact  if  we  proceed  as  above,  this  is  not  the  case.  To  see  this  consider 
the  sum 

and  writing  k  =(r-l)  we  have  , 

S(*V  j  2)  -  r^»  0  ^  rfk+of(V+o 

This  is  therefore  a  negative  binomial  series  and  hence 

n   ,     ,  ,-(A.)  VM 

Now  as  >-j>oo  ;  ?  ->  1  and  so  the  series  diverges.  However  if  we  exploit 
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the  terms  in  the  integral  outside  the  summation  sign,  we  can  circumvent 
this  problem.  We  consider  the  expression 

Now  since  * t  C"<vt]  We  have  ?■<  f0.TTe,]and  nence  the  series  does  not 
diverge.  Thus  replacing 

we  have  /  i  -, 

and  recalling  that  we  had  set  a  =  t1 

7      -     *~    v(T)       =     R^ 

which  fortunately  agree  with  the  results  obtained  earlier,  since  £"^f'M  -  \_ 
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